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Macroscopic parameters as well as precise information on the random force characterizing the 
Langevin type description of the nuclear fusion process around the Coulomb barrier are extracted 
from the microscopic dynamics of individual nucleons by exploiting the numerical simulation of the 
improved quantum molecular dynamics. It turns out that the dissipation dynamics of the relative 
motion between two fusing nuclei is caused by a non-Gaussian distribution of the random force. We 
find that the friction coefficient as well as the time correlation function of the random force takes 
particularly large values in a region a little bit inside of the Coulomb barrier. A clear non-Markovian 
effect is observed in the time correlation function of the random force. It is further shown that an 
emergent dynamics of the fusion process can be described by the generalized Langevin equation with 
memory effects by appropriately incorporating the microscopic information of individual nucleons 
through the random force and its time correlation function. 
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The fusion of two nuclei is one of the major non- 
equilibrium processes in low energy nuclear reactions 
where the fluctuation and dissipation play important 
roles. It is rather difficult to describe the fusion pro- 
cess without significant simplifications. Under various 
assumptions, several macroscopic transport models have 
been introduced to evaluate the formation of a compound 
nucleus in heavy-ion fusion reactions [1, 2]. However, the 
microscopic mechanism on how two colliding nuclei fuse, 
especially how the relevant kinetic energy dissipates into 
the intrinsic degrees of freedom (DoF), remains a subject 
requiring further research. 

On the other hand, it is becoming feasible to get var- 
ious information out of microscopic numerical simula- 
tions, like time-dependent Hartree-Fock (TDHF) theo- 
ries [3-7], the many-body correlation transport (MBCT) 
theory [8], the quantum molecular dynamics (QMD) [9], 
the antisymmetrized molecular dynamics [10], and the 
fermion molecular dynamics [11]. The TDHF theory is 
mainly based on the mean-field concept; in TDHF, fluc- 
tuations of collective variables are considerably underes- 
timated. Much effort has been made to give a beyond- 
mean-field description of fluctuations [12]. The n-body 
correlations are incorporated in the MBCT theory [8] 
which has only been used in very light systems [13]. 

The QMD is a microscopic dynamical n-body the- 
ory which was successfully used in intermediate-energy 
heavy-ion collisions (HIC) [9]. An improved QMD 
(ImQMD) has been developed in order to extend the ap- 
plication of QMD to low-energy HICs near the Coulomb 
barrier [14]. A series of improvements were made in the 
ImQMD,; in particular, by using the phase space occupa- 
tion constraint method [15], the fermionic properties of 


nucleons is remedied, which is important for low-energy 
collisions. Making full use of the microscopic informa- 
tion provided by ImQMD simulations, in this Letter, we 
try to understand how the macroscopic fusion dynamics 
emerges out of the microscopic one. 

We focus on a simplest case of symmetric fusion pro- 
cess with the impact parameter equal to zero. In this 
case, the system can be divided into the left- and right- 
half parts instead of a projectile and a target [16]. The 
relative motion between two centers of mass (CoM) of 
the left and right parts is chosen as the relevant DoF to 
be described by the Langevin equation. Our analysis is 
limited in a stage where the relative distance R is much 
larger than its width. 

The one-dimensional generalized Langevin equation 
with memory effects reads [17-19] 


where u(t) is the relative velocity between the two parts, 
oF (t) the random force felt by either part, u the reduced 
mass of the system, y(t—t’) the friction kernel, and V (R) 
the potential for the relative motion. 

In the ImQMD model [14], a trial wave function is re- 
stricted within a parameter space {r}, pj}, where r; and 
pj are mean values of position and momentum opera- 
tors of the jth nucleon which is expressed by a Gaussian 
wave packet. The time evolution of the trial wave func- 
tion under an effective potential is governed by the time- 
dependent variational principle [9-11]. An expectation 
value of the Hamiltonian is given by using an improved 
Skyrme potential energy density functional. In this Let- 
ter, we concentrate on head-on collisions of 9°Zr+?°Zr. 
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Ten thousand collision events were simulated. Each sim- 
ulation is started at R = Ro = 30 fm and with an incident 
energy E = 195 MeV. Numerical details can be found in 
Refs. [20]. 

The potential for the relative motion is defined as, 


V(R) = Etor(R) — Fiett(R) — Fright (R), (2) 
where Erot, Miett, and Erignt represent the energy of the 
system and those of the left and right parts, respectively; 
each of which consists of the kinetic energy, the nuclear 
and the Coulomb potential energies. The potential V (R) 
is shown in Figs. 1 and 3. The TDHF has also been used 
to extract microscopic interaction potentials between two 
nuclei [5, 21] which show similar features as those from 
the ImQMD simulations presented here and in Refs. [22]. 

The random force or the fluctuation of force in the ith 
event is defined as 


OF (x), = Fi(x) — (F(a2)), x=tor R, (3) 


where F;(x) = Sa fi (x) denotes the total force acting 
on the left (right) part of the system in the ith event, 
(F(x)) = 45; Fi(z) the mean value, and f} (x) the 
force on the jth nucleon in the left (right) part. Here 
A means the number of nucleons contained in the left 
(right) part and n denotes the total number of events. 
In Eq. (3) and hereafter, (Q) denotes an average of Q 
over all events. For low-energy collisions, the fluctua- 
tion mainly stems from the initialization of each event in 
which the position and the momentum of each particle 
are chosen randomly under certain conditions. With time 
this initial fluctuation propagates and is not smoothed 
out because in QMD a many-body rather than a mean 
field problem is solved [9]. 

Distributions of ôF (R) at several distances are shown 
in Fig. 1. The random force at R(t = 0) = Ro shows 
a Gaussian distribution with the full width at half max- 
imum (FWHM) r = 0.1 MeV/fm which could be un- 
derstood analytically as only the Coulomb field is felt 
by the particles. In a region far away from the barrier, 
e.g., R ~ 18 fm, ôF has a Gaussian distribution with 
Tr = 0.5 MeV/fm. From a certain distance, R ~ 13.5 
fm, there appears a non-Gaussian shape, as is observed 
in Fig. 1. According to the shape of the distribution 
of F(R), one may divide the whole process into three 
regions. Region 1 represents an approaching phase up 
to the touching point: The distribution has a Gaussian 
form with a rather narrow width. Region 2 is from the 
touching point to the barrier top: A non-Gaussian shape 
appears. Region 3 is from just inside the barrier top to 
the fusing phase: The distribution of F(R) has again 
a Gaussian shape with IT ~ 15 MeV/fm which is almost 
two orders of magnitude larger than that in Region 1. 

To make clear what happens in Region 2, we divide the 
distribution of ôF (R) into a symmetric Gaussian and an 
asymmetric tail parts as is shown in Fig 2. The width of 
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FIG. 1. (color online) Distributions of the random force 
F(R). Each inset shows the potential V (R) with the blue 


dot representing the position where the system locates. The 
contour plots display the nucleon density distribution of the 
system. 
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FIG. 2. (color online) Distribution of the random force ôF (R) 
at R = 13.5 fm which is divided into the symmetric Gaus- 
sian (dark blue) and asymmetric tail (light blue) parts. Two 
typical events are shown in the inset: The abscissa and the 
ordinate express relative position z of each nucleon and the 
force it feels in the z direction. 


the Gaussian part is of the same order of magnitude as 
that in Region 1. The detailed structure of the random 
force can be studied by examining the strength and direc- 
tion of the force felt by each nucleon. One typical event 
in the symmetric part is shown in Fig. 2(a): All nucleons 
are well divided into two separated groups expressing the 
projectile and the target, respectively. Moreover, each 
nucleon locating in the left side of each nucleus feels a 
force toward the right (positive value), and that in the 
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right side feels a force toward the left (negative value), so 
as to keep a stable mean-field. The resultant force made 
by all nucleons in each nucleus is almost zero. Namely, 
the intrinsic structure of two fusing nuclei is kept almost 
unchanged, so is the width of the random force. This sit- 
uation persists in events which belong to the symmetric 
Gaussian part in Region 2 and in all those in Region 1. 

A typical event in the asymmetric tail is shown in 
Fig. 2(b). Nucleons are roughly divided into two groups 
surrounded by solid lines. However, there appears a small 
third group within the dashed line. Since a few points in 
the negative (positive) force region express a set of nu- 
cleons which escape from the left (right) nucleus, and are 
being absorbed by the right (left) nucleus, a resultant 
force made by these nucleons gives a large right(left)- 
directed component to the random force. These trans- 
ferred nucleons move in an average potential formed by 
both the projectile and the target; they play a role to 
open a window. 

When the two nuclei come much closer, there occur 
more events which have more nucleons in the third group. 
Meanwhile, the other two groups, originating from the 
projectile and target, become closer to each other. Conse- 
quently, the asymmetric tail in the distribution of ôF (R) 
becomes larger. At the border between Regions 2 and 
3, it becomes very difficult to distinguish an event in 
the center part of the distribution from that in the tail 
part and all events are absorbed into a widely spreading 
Gaussian distribution. 

From above discussions, it is concluded that the main 
microscopic origin of the random force, i.e., a two orders 
of magnitude enhancement of the random force is gen- 
erated by individual nucleons in the third group. These 
nucleons also result in the abnormal behavior in the dis- 
tribution of 6F(R), i.e., the long tail in Region 2 and a 
much larger width in Region 3 compared to Region 1. 

Next let us extract information for the macroscopic 
dynamics out of microscopic simulations. Assuming 
that the work done by the friction force is completely 
converted into the intrinsic energy Eintr(R), one gets 
the friction coefficient yo(R) from the Rayleigh for- 
mula [6, 16, 19], 


(Frric( R)) 
Pr 4) 


with Frric(R) = dE intr (R)/dR, Eintr (R) = Erot (R) = 
Eeou(R), and Ecou(R) = P?/2u+V(R). P denotes 
the relative momentum between two CoMs and its mean 
value (P)p at a given R is defined as 


IDRO) 


where P;(t) and R;(t) are the momentum and coordinate 
of the i-th event at time t and the following correspon- 
dence is used: For each event i, a time t; is chosen in 


Yo(R) = 


Vice [Ri(ti)=R} > (5) 
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FIG. 3. (color online) The correlation function 
(6F(R)5F(R))*** [red dots, in (MeV/fm)?] and the friction 
coefficient yo(R) (blue squares, in 0.001 c/fm). The grey 
line shows the potential V(R). Pink diamonds represent 
(OF (R)OF(R))*™ calculated by eliminating events in the 
asymmetric tail. 


such a way that the relative distance takes a given value 
R, i.e., Ri(t;) = R. A R-dependent correlation function 
is defined as 


(OF (R)OF 


=- P 


(tile, |Ri(ti)=R}° 
(6) 


Figure 3 shows the correlation function (dF (R)dF(R)) 
and the friction coefficient yo( R) which play decisive roles 
in the macroscopic description of dissipation phenom- 
ena. As is seen from Fig. 3, (OF (R)OF(R)) and yo(R) 
have similar shapes and their peaks locate at similar R. 
The friction coefficient of the fusion process induced by 
a head-on collision extracted from TDHF calculations 
shows similar strong peak structure. As the incident en- 
ergy E increases, the shape of the curve yo(R) ~ R may 
change. When F is high enough, yo(R) increases gradu- 
ally with decreasing R [6, 16, 23). 

To explore more deeply the dynamical relation between 
the microscopic motion of individual nucleons and the 
macroscopic dissipative motion, in Fig. 4 we show the 
time correlation function of the random force o(R,7T) 
which is defined as 


= LY 6F y(t) oP iti 
ihe {=l 


T)linRrep=r} (7) 


In Fig. 4 one clearly finds the non-Markovian effect. 
Especially when R = 12 ~ 10 fm, it is important to take 
account of memory effects generated by the microscopic 
motion of nucleons when one tries to properly evaluates 
macroscopic effects of the dissipation. Starting from the 
generalized Langevin equation (1) with memory effects, 
one gets a generalized fluctuation-dissipation (GFD) re- 
lation (OF (t)dF (t — t')) = uTy(t — t') which properly 
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FIG. 4. (color online) Time correlation function o(r,7) (7). 
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FIG. 5. (color online) Effective temperature in the Markovian 
limit (a) and the Non-Markovian one (b). The blue line shows 


4/ Eintr/a with a = Atotal/4. 


takes account of the time correlation of the random force. 
There are many ways to define the temperature for com- 
pound nuclei (see, e.g., Ref. [24]). Here we define an 
effective temperature for colliding systems by applying 
the GFD relation, 


1 Co 
zo dro(R,T). (8) 


Thon—Markov R)= 
(7) HYo( 


The effective temperature Tt!  Warkoy aS Well as the 
one from the Markovian approximation TYP, are 
shown in Fig. 5. ./Hintr/a representing the temperature 
of a compound nucleus in the Fermi gas model is also 
shown as a reference. Although Tt) yparkoy and Thee ikov 
differ by one order of magnitude, they both show a peak 
around the range where the asymmetric tail appears in 
the distribution of ôF (R). These peaks are related to the 
fact that the relative motion for events in the asymmet- 
ric tail part of the F(R) distribution is strongly affected 
by a few transferred nucleons between two fusing nuclei, 
i.e., by those in the third group of Fig. 2(b). The macro- 
scopic dynamics of the relative motion described by the 
one-dimensional Langevin equation (1) is not appropri- 


ate in Region 2. In other words, the appearance of the 
non-Gaussian distributed random force indicates a neces- 
sity of introducing a new macroscopic DoF. Whether or 
not this new DoF may be related to the formation of a 
neck is an open question [1, 25]. 

After eliminating the events in the asymmetric tail in 
the distribution of ôF (R), one gets effective temperatures 
TOM Markov 224 Tytarkoy Which are depicted in Fig. 5. 
The correlation function (F (R)ôF(R)) after eliminating 


ge ei. DA sym 
the asymmetric tail is also shown in Fig. 3. T on- Markov 


shows a consistent feature with y/Eintr/a in Region 3. 
While TÈ koy 18 by an order of magnitude smaller than 
TO Markov: That is, the amount of energy dissipated 
from the relative motion into the intrinsic DoFs could 
be more properly described by the generalized Langevin 
equation with memory effects. 

When the incident energy E is far above the Coulomb 
barrier, the non-Gaussian fluctuation and the non- 
Markovian effect become less pronounced [23]. It will be 
interesting to study the dependence of the non-Gaussian 
fluctuation and the non-Markovian effect on FE as well 
as the impact parameter and the reaction system. The 
spin-orbit coupling is important to properly reproduce 
the dissipation in heavy-ion fusion reactions [26]; e.g., 
the so-called “fusion-window” problem was solved in the 
first quantitative TDHF calculations with the inclusion 
of the spin-orbit interaction [26]. One may expect more 
dissipations if the spin-orbit coupling effects are included 
in the ImQMD simulations. 

In summary, we have discussed the generalized 
Langevin dynamics with memory effects by using both 
the macroscopic and microscopic information extracted 
from ImQMD simulations for the fusion process around 
the Coulomb barrier. It is found that the dissipation 
dynamics of the relative motion between two fusing nu- 
clei is associated with non-Gaussian distributions of the 
random force. In addition to the macroscopic informa- 
tion like the friction coefficient and the potential for the 
relative motion, the microscopic information of the ran- 
dom force as well as of its time correlation function and a 
proper treatment of the non-Markovian (memory) effect 
in the Langevin dynamics are decisive for the dynamics 
of emergence in the nuclear dissipative fusion motion. 
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